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ABSTRACT 

We perform a detailed and systematic investigation of the possible impacts of dark stars upon the reionization 
history of the Universe, and its signatures in the cosmic microwave background (CMB). We compute hydrogen 
reionization histories, CMB optical depths and anisotropy power spectra for a range of stellar populations 
including dark stars. If dark stars capture large amounts of dark matter via nuclear scattering, reionization can 
be substantially delayed, leading to decreases in the integrated optical depth to last scattering and large-scale 
power in the EE polarization power spectrum. Using the integrated optical depth observed by WMAP7, in our 
canonical reionization model we rule out the section of parameter space where dark stars with high scattering- 
induced capture rates tie up > 90% of all the first star-forming baryons, and live for > 250 Myr When nuclear 
scattering delivers only moderate amounts of dark matter, reionization can instead be sped up slightly, modestly 
increasing the CMB optical depth. If dark stars do not obtain any dark matter via nuclear scattering, effects upon 
reionization and the CMB are negligible. The effects of dark stars upon reionization and its CMB markers can 
be largely mimicked or compensated for by changes in the existing parameters of reionization models, making 
dark stars difficult to disentangle from astrophysical uncertainties, but also widening the range of standard 
parameters in reionization models that can be made consistent with observations. 

Subject headings: dark ages, reionization, first stars - dark matter - cosmic background radiation - stars: 
Population III 



1. EVITRODUCTION 

At a redshift of z ^ 1 100, the bulk of electrons and protons 
in the Universe recombined to form neutral hydrogen. This 
cleared the way for thermal photons to free-stream away from 
the resulting surface of last scattering, forming what we now 
see as the cosmic microwave background (CMB; see e.g. Hu 
& Dodelson 2002; Samtleben et al. 2007). A prolonged pe- 
riod of darkness ensued, until the first sources of hard ionizing 
radiation forming within galaxies appeared at redshifts z < 30 
(Gnedin 2000; Ciardi et al. 2000; Bromm et al. 2001; Wyithe 
& Loeb 2003; Schaerer 2002; Venkatesan et al. 2003; Benson 
et al. 2006; Loeb 2009). These sources are thought to have 
reionized the neutral gas in the intergalactic medium (IGM). 
This process was completed by z ^ 7 (Fan et al. 2006; Daw- 
son et al. 2007), with the IGM subsequently remaining in the 
fully-ionized state we see it in today. 

The process of reionization is expected to have left its im- 
print upon a number of cosmological observables. The CMB 
is sensitive to the total optical depth to the surface of last scat- 
tering; any free electrons between us and the surface will have 
scattered CMB photons before they could reach us, modify- 
ing the anisotropy power spectra we observe today. Similarly, 
the changing distribution of neutral hydrogen during reioniza- 
tion can be mapped using its ground-state hyperfine transition 
(Furlanetto et al. 2006; Chen & Miralda-Escude 2008), corre- 
sponding to a rest-frame wavelength of 21 cm. The signal is 
expected to be weak because this is a forbidden line, so the 
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transition probability is small. A multitude of upcoming ex- 
periments hope to detect it nonetheless (Morales & Wyithe 
2010); these include the Low-Frequency (LOFAR; Harker 
et al. 2010), Murchison Wide-Field (MWA; Mitchell et al. 
2010) and Square Kilometer Arrays (SKA; Carilli 2008). 

The first ionizing sources are generally thought to be 
the first stars (see, e.g., Gnedin 2000; Bromm et al. 2001; 
Venkatesan et al. 2003; Tumlinson et al. 2004; Wise & Abel 
2008), referred to as Pop III. The terms Pop III and Pop II 
are typically interpreted as broad distinctions based on com- 
position - respectively, they represent stars that are metal-free 
and metal-poor. Pop III stars consist entirely of primordial 
hydrogen and helium synthesized in the Big Bang, so are 
likely to have a mass function weighted towards higher stel- 
lar masses than that of Pop II, due to the absence of metal 
lines, which allow efficient gas cooling and cloud fragmenta- 
tion in metal-enriched galaxies at later epochs (Bromm et al. 
2001; Tumlinson et al. 2003). Within the Pop III category, the 
terms Pop III. 1 and Pop III. 2 have recently arisen in the lit- 
erature. Typically, Pop III. 1 connotes the very first stars (i.e., 
first-generation metal-free stars) that have masses exceeding 
-100- 300 Mq, and that form in -10^ M© dark matter (DM) 
minihalos. Pop III. 2 stars form in the wake of the radiative and 
chemodynamic feedback of Pop III. 1 supernovae. Pop III. 2 
are therefore "second-generation" stars with metal-free com- 
position, and are thought to have lower stellar masses on av- 
erage (—10- IOOM0) than Pop III.l stars (Tumlinson et al. 
2004; McKee & Tan 2008; Ohkubo et al. 2009; see, how- 
ever, Clark et al. 201 1 for arguments on Pop III. 1 stars having 
lower stellar mases than Pop III. 2 stars). Although Pop III.l 
stars were initially thought to form in isolation, producing a 
single very massive star per minihalo (Abel et al. 2002), this 
paradigm has lately given way to one where they form mostly 
in pairs, or systems of even higher multiplicity (Krumholz 
et al. 2009; Turk et al. 2009; Stacy et al. 2010). 

Recent observations of the galaxy luminosity function (LF) 
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at high redshifts, z = 4-8 (Bouwens et al. 2010, 2011) indi- 
cate that the earliest, faintest galaxy halos make a substantial 
contribution to reionization, in good agreement with the hy- 
pothesis that the first stars are indeed the first ionizing sources. 
Connecting stellar population models to UV observations of 
the faint-end galaxy LF is fraught with systematic uncertain- 
ties, however, as it requires extrapolation of the observations 
to even higher redshifts (z > 8), and making a number of as- 
sumptions with regard to the parameterization of the under- 
lying reionization models. For these reasons, and given the 
relative insensitivity of present-day 21 cm data to reionization 
physics, we focus on constraints from the CMB in this paper. 

It has been shown (Spolyar et al. 2008; Natarajan et al. 
2009) that dark matter could have had a more direct impact 
on Pop III.l star formation than simply facilitating the ini- 
tial baryonic collapse. If DM consists of a new particle that 
is present in equal numbers to its antiparticle, or if it is in- 
deed its own antiparticle, it will self-annihilate to produce 
Standard Model (SM) particles such as quarks, photons and 
electrons. This is the case for e.g. weakly-interacting massive 
particles (WIMPs), arguably the most widely-studied and nat- 
ural solution to the dark matter problem (see e.g. Jungman 
et al. 1996; Bergstrom 2000; Bertone et al. 2005; Bertone 
2010). As baryons cool and contract during star formation, 
they steepen the gravitational potential within the minihalo, 
drawing even more dark matter into its center (Freese et al. 
2009). The resultant spike in the DM annihilation rate can 
inject an appreciable amount of energy into the collapsing 
cloud, halting or delaying star formation (Mapelli et al. 2006; 
Stasielak et al. 2007; Ripamonti et al. 2007, 2010) and re- 
sulting in a cool, partially-collapsed object known as a "dark 
star" (Spolyar et al. 2008; Natarajan et al. 2009). Ongoing 
annihilation of DM particles in the core of a star can have 
substantial impacts upon its structure and evolution (Salati & 
Silk 1989; Fairbairn et al. 2008; Scott et al. 2008; locco 2008; 
locco et al. 2008; Freese et al. 2008a; Scott et al. 2009; Spol- 
yar et al. 2009). 

As dark stars have significantly different structures and evo- 
lutionary histories to "normal" Pop III stars, their ionizing 
photon outputs differ substantially as well (Yoon et al. 2008), 
leading to a potentially distinct impact upon the process of 
reionization (Schleicher et al. 2009). The purpose of this pa- 
per is to systematically investigate the effects of dark stars 
upon the reionization history of the Universe. We begin by 
discussing dark star formation and evolution in Sec. 2, and 
stellar population models including dark stars in Sec. 3. We 
give an overview of our reionization models and calculations 
in Sec. 4. In Sec. 5 we give the resulting alternative reion- 
ization histories for Universes containing dark stellar popula- 
tions. In Sec. 6, we show how such IGM ionization histories 
would impact the measured optical depth to electron scatter- 
ing in the CMB, drawing limits on dark star populations using 
the integrated optical depth measured by WMAP7, and mak- 
ing predictions for the corresponding constraining power of 
Planck. We also show the impact of dark stars upon the EE 
polarization anisotropy spectrum of the CMB, and discuss its 
potential use for constraining dark star populations. We con- 
clude in Sec. 7. 

2. FORMATION AND EVOLUTION OF DARK STARS 

Assuming DM self-annihilates, the formation of a dark star 
relies on some efficient means for bringing the DM into the 
center of a star The two processes which may provide this 
means are gravitational contraction and nuclear scattering. 



The first is simply an effect of the changing gravitational 
potential during the collapse of a baryonic gas cloud (Freese 
et al. 2008b). As baryons cool and collapse onto the central 
overdensity, dissipating energy by radiative emission and an- 
gular momentum by hydrodynamic and magnetic interactions 
(which may be enhanced during the formation of the first stars 
by dynamo effects; Schleicher et al. 2010; Sur et al. 2010), the 
gravitational potential in the core of the halo steepens. In turn, 
the steepening potential draws DM into the center of the cloud 
(despite its inability to actually dissipate energy), resulting in 
a strongly-peaked DM density distribution. This is referred 
to simply as 'gravitational contraction,' and is more general 
than the well-known case of adiabatic contraction because it 
does not strictly require that the gravitational potential change 
more slowly than the orbital timescale of individual particles. 
In the canonical scenario, where a single Pop III. 1 star forms 
at the center of the very first halos, stars begin their lives al- 
ready in possession of a large reservoir of DM. 

The second way for DM to end up in a stellar core relies on 
it possessing a weak-scale scattering cross-section with nu- 
cleons. Such an interaction is characteristic of WIMP dark 
matter. Assuming such a cross-section, DM particles passing 
through a star can lose energy through collisions with stellar 
nuclei, becoming gravitationally bound to the star (Steigman 
et al. 1978; Krauss et al. 1985; Press & Spergel 1985; Gould 
1987). This leads to repeat scattering events, eventually re- 
moving enough energy that the particle ends up in the stellar 
core. 

Regardless of the path DM follows into a star, the effects 
are essentially the same. DM annihilation in the core provides 
an additional energy source alongside nuclear fusion, causing 
the core to expand and cool. This occurs due to the negative 
specific heat of a self-gravitating body, and the fact that the 
DM annihilation rate is decoupled from the nuclear core den- 
sity. The core expansion leads to a larger, cooler, typically 
strongly convective stellar object (Spolyar et al. 2008; Scott 
et al. 2009; Spolyar et al. 2009; Casanellas & Lopes 2009). 
In the case of gravitationally-contracted DM, the collapse of 
the forming star is slowed, effectively extending the proto- 
stellar phase. The slowdown allows gas accretion to continue 
longer than it otherwise would, as the onset of radiative feed- 
back is delayed. The resulting object therefore grows more 
than in the absence of DM (Spolyar et al. 2008; Umeda et al. 
2009), leading to masses of order --800-lOOOMo. It has been 
suggested that supermassive objects might even be possible 
(Freese et al. 2010), though this is strongly constrained by 
existing data (Zackrisson et al. 2010b). 

The degree to which a star affected by DM annihilation re- 
sembles either a fully-fledged dark star (during the extended 
protostellar phase) or a main-sequence object depends upon 
the rate of annihilation in its core. Higher rates of annihilation 
are required to support larger, more diffuse, protostellar-like 
structures against further collapse. The evolution of a dark 
star therefore depends strongly on the rate at which dark mat- 
ter is delivered to the stellar core, and how that rate changes 
over time. Typically, gravitationally-contracted dark matter 
will be exhausted in a period of ^0.4 Myr (Spolyar et al. 
2009).^ Without replenishment via DM capture due to nu- 

^ We note however that some uncertainty remains over this value, with 
much shorter timescales suggested by locco et al. (2008). Ripamonti et al. 
(201 0) even show that at least in the early stages of the collapse, dark matter 
annihilation might in fact help the gas to cool and contract rather than hinder 
it, by enhancing the formation of H2 molecules; this effect has yet to be taken 
into account in most dark star modeling. 
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clear scattering, dark stars then contract, heat up and move 
on to the main sequence (MS) to live their Uves as extremely 
massive Pop III.l stars. 

In the simplest scenario, where the DM halo is homoge- 
neous and spherically symmetric, DM capture by nuclear 
scattering has been shown not to substantially extend the life- 
time of dark stars (Sivertsson & Gondolo 2011). The im- 
pact of more realistic halo distributions upon the DM star- 
crossing rate, and therefore the capture rate, remains to be 
understood however (see e.g. Freese et al. 2010 for the sug- 
gestion that particle orbits may even be strongly centrophilic). 
In principle, if capture rates due to nuclear scattering are high 
enough, dark stars may exist as cool, diffuse objects for up to 
^500 Myr (see e.g. Zackrisson et al. 2010a, for a discussion 
of possible dark star lifetimes). 

Similarly, the impact of Pop III. 1 stars forming as binaries 
or higher-multiplicity systems is not yet well understood. At 
some level the displacement of the collapsing baryonic core(s) 
from the central DM spike may indeed prevent dark star for- 
mation all together. Alternatively, the fragmentation process 
may introduce sufficient structure to the phase-space distribu- 
tion of the DM halo that the findings of Sivertsson & Gon- 
dolo (2011) are circumvented, and dark star lifetimes in fact 
become longer. In this case however, the resulting dark stars 
might be of substantially lower mass than those formed in iso- 
lation, due to reduced accretion onto the central object (Peters 
et al. 2010), shown to be significant for first star formation 
(Clark et al. 201 1). Further review of the structure and evolu- 
tion of dark stars and the mechanisms for fueling them can be 
found in Scott (2011). 

Given the substantial theoretical uncertainty in predicting 
the longevity of dark stars, for the purposes of this paper we 
consider their lifetimes to be a free parameter, to be con- 
strained by CMB data or other observations of reionization. 
Although we choose a halo mass threshold at z = 20 that is de- 
signed to allow the formation of a single dark star in the very 
smallest halos (see Sec. 4), we do not make any strong as- 
sumptions as to the number of dark stars forming in each halo 
in general, as the reionization formalism described in Sec. 4 
is non-specific as to the multiplicity of star formation. 

3. DARK AND SEMI-DARK STELLAR POPULATIONS 
3.1. Population models 

Modeling a stellar population containing dark stars requires 
careful consideration of the possible evolutionary histories of 
the dark component. We define different structural forms for 
dark stars, depending upon the degree to which they are dom- 
inated by dark matter: the dark star proper (DSP) and the dark 
star near the main sequence (DSNMS). The DSP structure oc- 
curs when dark matter annihilation contributes a substantial 
fraction of the star's energy budget, leading to the large, dif- 
fuse, cool objects discussed by Spolyar et al. (2008). These 
objects lie far to the right of the HR diagram, in the region 
populated by protostars on their way to the main sequence. 
The DSNMS structure occurs when dark matter contributes 
only a small amount of the star's total energy budget, and 
much more closely resembles a main sequence star. Corre- 
spondingly, they lie only slightly to the right of the standard 
Pop III main sequence on the HR diagram. 

Of course, the two structures are not entirely distinct; a con- 
tinuum of stellar structures is possible, parameterized by the 
amount of energy produced by dark matter annihilation in the 
stellar core (as discussed in Sec. 2). Exactly what structure 



a dark star exhibits depends upon its age (with the DSNMS 
stage always following the DSP stage) and the rate at which 
dark matter is captured and converted into heat. Sustained 
high rates of capture therefore effectively increase the dura- 
tion of the DSP phase, but not the DSNMS phase, as the star is 
supported by DM annihilation in a cool, diffuse configuration 
until the DM runs out. Correspondingly, sustained moderate 
capture rates increase the duration of the DSNMS phase, but 
not the DSP phase, as the star is supported by DM annihilation 
and nuclear burning, in a configuration that is only slightly 
cooler than the equivalent main sequence structure (see e.g. 
Scott et al. 2009). 

We take a phenomenological definition for the exact de- 
marcation between the DSP and DSNMS phases. The DSP 
phase includes all structures where capture rates are suffi- 
ciently large to keep the star too cool to contribute to reioniza- 
tion {Q ^ 0; cooler stars are redder, so more of their luminos- 
ity is output at wavelengths longward of the neutral hydrogen 
ionization threshold). Structures in the DSNMS phase, on the 
other hand, have small enough WIMP capture rates to be suf- 
ficiently hot to make some contribution to reionization. Be- 
cause the ionizing photon flux for a given stellar mass falls 
off very abruptly as the relative contribution of dark matter 
annihilation to the star's total energy budget is increased, to 
a first approximation stars in this phase can be modeled as 
main-sequence objects. 

We define three distinct capture scenarios: no capture (NC), 
meager capture (MC) and extreme capture (EC). NC is the 
canonical scenario discussed and simulated using a simple 
varying-index polytropic model by Spolyar et al. (2009): dark 
stars form by gravitational collapse, and grow to exhibit larger 
masses than typical Pop III.l stars because annihilation of 
gravitationally-contracted dark matter in their cores inhibits 
the collapse. After the initial (gravitationally-contracted) pop- 
ulation of dark matter is exhausted, the star finishes contract- 
ing and makes its way to the main sequence, where it lives like 
any other ^^00 Mq star. In this case, we have a DSP phase 
with a duration approximately equal to the ages of the stars of 
Spolyar et al. (2009) upon reaching the MS (fosp ^ 0.4 Myr 
for a DM mass of m^ = 100 GeV), followed by a DSNMS 
phase of duration equal to the standard MS lifetime of an 
--800 Mq MS Pop III.l stai- (fosNMS ^ 1 Myr; e.g. Schaerer 
2002). 

In the MC scenario, we have enough dark matter capture 
to extend the lifetime fosNMS of the DSNMS phase, whereas 
in the EC scenario we have enough dark matter capture to 
instead extend fosp, the lifetime of the DSP phase. In these 
two scenarios, we therefore take the lifetimes of the respective 
phases as free parameters, and keep the lifetime of the other 
phase fixed at its canonical value in the NC scenario. 

In general, the maximum duration of any phase of the life- 
time of a star powered by dark matter capture (either wholly 
or partially) is limited by the total core hydrogen-burning life- 
time of the star, and the self-annihilation time of the dark mat- 
ter halo from which the star captures its dark matter. For stars 
in the most dark-matter dominated parts of the DSP phase, the 
low core temperature and very high capture rate mean that the 
self-annihilation time is the relevant limit, and the hydrogen- 
burning time plays little role. For stars in the DSNMS phase, 
the opposite is true. These limits were calculated and dis- 
cussed in detail by Zackrisson et al. (2010a). Here we approx- 
imately adopt the values fmax of Zackrisson et al. as the upper 
limits for the durations of the two phases. The lower limits are 
given by the NC case, which can be seen as a limiting case of 
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TABLE 1 

Scenarios, parameters, redshifts of H i reionization 

AND integrated OPTICAL DEPTHS TO LAST SCATTERING 
FOR STELLAR POPULATIONS WITH DARK STARS. OPTICAL 
DEPTHS INCLUDE CONTRIBUTIONS FROM HE II, HE III AND 
RESIDUAL ELECTRON FRACTION AFTER RECOMBINATION. 
See TEXT FOR DETAILS. 



Scenario 


fDSP 
(Myr) 


fDSNMS 
(Myr) 


/ds 


^reion 


Te 


Pop II+III only 











11.050 


0.1084 


NC 


0.4 


1 


0.01 


11.050 


0.1084 


(no capture) 






1 


11.050 


0.1084 


MC 


0.4 


3 


0.01 


11.050 


0.1084 


(meager capture) 






0.1 


11.061 


0.1085 








0.6 


11.092 


0.1089 








1 


11.122 


0.1093 






6 


0.01 


11.050 


0.1084 








0.1 


11.071 


0.1086 








0.3 


11.102 


0.1091 








0.6 


11.143 


0.1097 








0.8 


11.174 


0.1101 








1 


11.205 


0.1105 


EC 


5 


1 


0.01 


11.050 


0.1084 


(extreme capture) 






0.1 


11.030 


0.1080 








0.3 


10.978 


0.1072 








0.6 


10.885 


0.1060 








0.8 


10.822 


0.1051 








1 


10.760 


0.1042 




15 


1 


0.01 


11.050 


0.1083 








0.1 


10.999 


0.1076 








0.3 


10.885 


0.1059 








0.6 


10.697 


0.1031 








0.8 


10.560 


0.1011 








0.9 


10.486 


0.0999 








1 


10.423 


0.0988 




50 


1 


0.01 


11.040 


0.1083 








0.1 


10.968 


0.1072 








0.3 


10.780 


0.1046 








0.6 


10.476 


0.1000 








0.8 


10.231 


0.0962 








0.9 


10.103 


0.0941 








1 


9.973 


0.0917 




150 


1 


0.01 


11.040 


0.1083 








0.1 


10.936 


0.1069 








0.3 


10.644 


0.1033 








0.6 


10.060 


0.0957 








0.8 


9.514 


0.0878 








0.9 


9.192 


0.0824 








1 


8.852 


0.0755 




500 


1 


0.01 


11.040 


0.1083 








0.1 


10.926 


0.1069 








0.3 


10.623 


0.1032 








0.4 


10.433 


0.1010 








0.5 


10.210 


0.0983 








0.6 


9.941 


0.0951 








0.7 


9.592 


0.0910 








0.8 


9.079 


0.0852 








0.85 


8.726 


0.0811 








0.9 


8.215 


0.0755 








1 


6.282 


0.0446 



both the MC and EC scenarios (although physically, it is of 
course only continuously linked to the MC case). This gives 
0.4 Myr < tusp < 500 Myr and 1 Myr < fosNMS ^ 6 Myr. 

Finally, we also have the dark star mass fraction /ds as an 
additional free parameter in all three scenarios (NC, MC and 
EC). This describes the fraction of the star-forming baryonic 
mass that initially goes into dark stars rather than normal Pop 
III stars. We begin with an initial population consisting of a 
fraction /ds of dark stars in the DSP phase and a complemen- 
tary fraction (1 -/ds) of normal Pop III stars. As the normal 
Pop in stars finish their starburst (i.e. after fpopiii, which we 



set to 10 Myr), they are replaced with Pop II stars. As the 
dark stars transition from DSP to DSNMS and then eventu- 
ally die, they are replaced with either normal Pop III stars if 
the death occurs before the end of the original Pop III star- 
burst, i.e. fDSP+fDSNMS < fpopiii = lOMyr, or directly with Pop 
II stars if the death occurs after the cessation of the original 
Pop III starburst, i.e. Jdsp + fDSNMS > fpopiii = 10 Myr. In the 
case where the dark stars die and are replaced by a new popu- 
lation of Pop III stars, these Pop III stars later also die and are 
replaced by Pop II stars at f = fpopin = 10 Myr after the begin- 
ning of star formation, just like their counterparts in the com- 
plementary fraction (1 -/ds) of normal Pop III stars present 
from the beginning of the calculation. 

For ease of reference, the full set of scenarios and parame- 
ters with which we compute reionization histories is given in 
Table 1. 

For some combinations of parameters (large /ds and long 
fDSp). the astute reader will have realized that our population 
models contain either very few or no normal Pop III stars. 
Constructing a consistent picture of the chemical evolution of 
the Universe in these cases becomes somewhat more problem- 
atic than in the standard situation, where supernovae produced 
by the deaths of the original Pop III stars provide the neces- 
sary chemical enrichment of gas to facilitate the formation of 
Pop II stars. Stars as heavy as the 800 M© dark stars we con- 
sider here are typically expected to collapse directly to black 
holes (e.g. Umeda et al. 2009), producing very few metals. 
Because only stars in the mass range up to 260 Mq produce 
metals (Heger & Woosley 2002; Venkatesan & Truran 2003), 
this constrains the mass range of a primordial stellar popula- 
tion that must necessarily seed the conditions for Pop II star 
formation (and this is used to justify the Pop III mass func- 
tion we consider below). Moreover, the yields of heavy el- 
ements from supernovae from the first stars as well as the 
true distribution of metallicities in the very first Pop II stars 
are currently not well-constrained, so it is conceivable that a 
small number of highly-efficient Pop III stars could provide 
the bulk of chemical enrichment necessary to allow the tran- 
sition to Pop II (and such a Pop II may have begun somewhat 
later in some locations, and with a somewhat smaller metal- 
licity than is considered typical). In any case, stars as massive 
as 800 M0 operate very close to the Eddington luminosity, so 
are expected to exhibit very strong stellar winds and experi- 
ence numerous mass-loss events. Although mass loss from 
metal-poor or metal-free stars is expected to be substantially 
reduced in comparison to that from their metal-rich cousins, 
depending upon the rotational and convective properties of the 
first stars, the material blown off from such objects may well 
be sufficiently processed to also contribute to the chemical 
enrichment of the Universe. 

3.2. Ionizing photon fluxes 

At each timestep of our reionization calculation, we cal- 
culate the weighted-average, mass-normalized, hydrogen- 
ionizing photon output of the combined stellar population as 



etot(0 = /DSeDs(0 + (l -/Ds)enormal(0- 



(1) 



Here 2ds(0 refers to the ionizing photon output per unit mass 
of stars in the population originally consisting of dark stars. 
Depending upon the time t in question, and the values of the 
lifetime parameters fDSP and Tdsnivis, this may be equal to ei- 
ther Qdsp, Gdsnms, Gpopii or Qpopiii- Gnormai(0 rcfcrs to the 
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population originally consisting of normal Pop III stars, and 
is equal to Qpopiii for t < lOMyr, and Qpopii for t > lOMyr. 

Explicitly, if we designate to as the time elapsed since the 
onset of star formation (DS and/or Pop III), and define the 
specific averaged Qu factors 

2tot.DSP+P3 = (l-/Ds)2popni (2) 

Ctot,DSP+P2 = (1-/ds)2popII (3) 

Gtot,DSNMS+P3 = /dsGdSNMS+(1~/ds)2popIII (4) 

2totDSNMS+P2 = /ds2dSNMS + (1-/ds)2popII, (5) 

we have three possible scenarios, depending upon the relative 
values of the lifetime parameters fosp, ^dsnms and fpopiii: 

• foSP + foSNMS < f Pop III 

1- 2h = Gtot,DSP+P3, duration fDSP 

2- Qh = 2tot,DSNMS+P3, duration fnsNMS 

3- Gh = Gpopiii, duration rpopiii-(fDSNMS +fDSp) 
4. ge = fipopii, duration fo-fpopiii 

• fosp < f Pop III and fosp + ^DSNivis > fpopiii 

1- 6h = 2tot,DSP+P3, duration fDSP 

2. Qh = 2tot,DSNMS+P3, duration fpopiii-fosp 

3- 2h = 2tot,DSNMS+P2, duration fnsNMS - (fpopiii - 
^Dsp) 

4. 6h = fipopii, duration fo-(fDSP + fDSNMs) 

• foSP > f Pop III 

1- 2h = 2tot,DSP+P3, duration fpopiu 
2. 2h = Gtot,DSP+P2, duration fosp-fpopiii 
3- 2h = Gtot,DSNMS+P2, duration fDSNMS 
4. Qh = Qpopii, duration fo-(fDSP + fDSNMs) 

For all scenarios, and in Eqs. 2 and 3 above, we assume 
Gdsp = during the DSP phase. 

We calculate ionizing photon fluxes during the DSNMS 
phase using model dark star atmospheres computed with 
TLUSTY (Hubeny & Lanz 1995), as described in Zackris- 
son et al. (2010a). Hydrogen-ionizing photon fluxes for some 
example dark star models computed by Spolyar et al. (2009) 
are shown in Table 2. Here we display Qu values for three dif- 
ferent stellar models, computed assuming three different DM 
masses. We show snapshots of Qh at different times in the re- 
spective models' evolution: larger stellar masses correspond 
to later models, as dark stars gradually accrete more matter. 
We do not show Qu values for earlier times in the simulations 
(corresponding to the main part of the DSP phase), as the ion- 
izing fluxes of the earliest (lowest-mass) snapshots shown in 
Table 2 are already too low to be significant for reionization. 
As is to be expected, ionizing photon fluxes increase with time 
and stellar mass as the dark stars become hotter, more com- 
pact and luminous as they move from DSP to DSNMS, and 
finally, to the zero-age main sequence. 

For a fixed initial DM density, larger DM masses lead to 
slightly decreased energy production in the stellar core (due 
to the decreased number density implied by a constant mass 
density). This can be seen in the lower final mass of the 
m^ = 1 TeV model in Table 2, where DM annihilation has ex- 
tended the accretion phase during the DSP less than it would 



TABLE 2 

Ionizing photon fluxes, with Qh values high enough to 

contribute to reionization (cf. table 1 in tumlinson 

et al. 2004) marked in bold. 



DM Mass 


Age 


Stellar 


Qh 


Qh 




(yr) 


Mass (Mq) 


(s-') 


(s-i MqI) 


IGeV 


3.1 X 10^ 


756 


1.05 xlO'*'' 


1.38x10"' 


IGeV 


3.3 X 10^ 


793 


4.26 X lO'*' 


5.37x1046 


IGeV 


4.6 X 10^ 


824 


5.37x1050 


6.52x10'" 


100 GeV 


3.0 X 10^ 


716 


6.28 xlO'*^ 


8.78x10''** 


100 GeV 


3.9 xlO^ 


779t 


4.97x1050 


6.38x10*" 


100 GeV 


4.1 X 10^ 


787* 


5.20X lO'O 


6.61 X 10"' 


10 TeV 


0.9 X 10^ 


327 


3.40X 10« 


1.04x10"' 


10 TeV 


2.7 X 10^ 


553 


3.27x10^0 


5.92x10*" 


Normal Pop 


11 


1-100 




7.76 xl0"6 


Normal Pop 


111 


10-140 




4.30x10*" 



' This case corresponds to the NC scenario considered in this work, 
and does not include any capture by nuclear scattering. All other en- 
tries in this table, including the model used for the MC and EC scenar- 
ios, correspond to models where nuclear scattering provides a similar 
amount of power to DM obtained by gravitational contraction (see 
Spolyar et al. 2009, for details). 

* This case is used in calculations for the EC and MC scenarios con- 
sidered in this work. 

if the DM mass were smaller (see Spolyar et al. 2009, for 
more details), resulting in a reduced Q^ during the DSNMS 
phase. In general however, the mass of the DM particle has 
only a weak impact upon the phenomenology of dark stars. 
For the remainder of this paper, we focus on the example case 
mx = 100GeV. 

For the NC scenario, we use the final, 779 Mq, m^ = 
100 GeV stellar model of Spolyar et al. (2009), which was 
computed without including any capture of WIMPs via nu- 
clear scattering. When calculating ionizing fluxes during the 
DSNMS phase in the MC and EC scenarios, we instead use 
the corresponding 787 M©, m^ = 100 GeV model, which in- 
cluded a small amount of DM capture by nuclear scattering 
(see Spolyar et al. for details). In practice, there is very little 
difference between the models of Spolyar et al. (2009) with 
and without the very small amount of capture they included. 

The ionizing photon fluxes of our canonical Pop II and III 
populations are also given in Table 2. 

4. REIONIZATION CALCULATIONS 

We use the semianalytic reionization model in Venkatesan 
et al. (2003) for a ACDM cosmology. The growth of ionized 
regions is tracked by a Press-Schechter formalism in combi- 
nation with numerical solutions for the growth of individual 
ionization fronts. We take our cosmological parameter set 
from the latest WMAP7 results (Larson et al. 2011).** We 
assume that the fraction of baryons forming stars in each halo 
is /i = 0.05, and that the escape fraction of H II ionizing ra- 
diation from halos is /esc = 0.1. Theoretical calculations, both 
analytical and from numerical simulations, as well as observa- 
tions of star-forming galaxies in the local and high-z Universe, 
indicate that /^ = 0.01-0.1, and f^sc = 0.01-0.2 (Venkatesan 
et al. 2003, and references therein). 

We allow star formation in all halos of virial temperature 
> 375 K (rather than 10"^ K as in Venkatesan et al. 2003, but 

^ We used the WMAP7 cosmological parameter set at: http://lambda.gs- 
fc.nasa.gov/product/map/current/params/lcdm_sz_lens_wmap7.cfm 
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similar to lO-' K in Tumlinson et al. 2004) starting at z = 20. 
This ensures that in our cases where /ds = 1, there are suffi- 
cient baryons in a 10^ M0 dark matter halo at z ^ 20 for our 
adopted star formation efficiencies to form a single 800 Mq 
dark star. This relies on the assumption that halos with virial 
temperatures < 10"* K can cool effectively using, e.g., H2 or 
other pathways (Abel et al. 2002; Haiman & Holder 2003). 
"^ Increasing the virial temperature threshold (i.e. the min- 
imum mass scale of collapsing halos) to values higher than 
our choice here will delay reionization correspondingly. As 
the most interesting constraints we produce are based on dark- 
star induced delays to reionization, this choice is a conserva- 
tive one. 

We assume Pop III (metal-free) stars can form starting at 
z =20, in a Salpeter-slope initial mass function (IMF) span- 
ning the range 10-140 Mq. This Pop III IMF'" is one lacking 
in low -mass stars, and was motivated by ionization constraints 
and observations of abundance trends in metal-poor Galactic 
halo stars (Tumlinson et al. 2004). It has a mass-normalized 
ionizing photon flux about a factor of 3 smaller than the ion- 
izing flux of an IMF containing very massive stars (lO^-lO^* 
M0, Bromm et al. 2001). We assume that the Pop III phase 
lasts for lOMyr in each halo, in agreement with the duration 
of metal-free star formation calculated from numerical sim- 
ulations of halo self-enrichment (Wada & Venkatesan 2003), 
and of interhalo enrichment (Bromm et al. 2003; Tumlinson 
et al. 2004). Subsequent to lOMyr after the onset of star for- 
mation in each halo, the ionizing spectrum is switched to a 
representative example of Pop II stars in a Salpeter-slope IMF 
in the mass range 1-100 Mq with metallicity Z = 0.001 (Lei- 
thereretal. 1995). 

It is possible that Pop II star formation need not always 
follow that of Pop III; both populations could form simulta- 
neously in separate gas clouds within individual halos, due 
to inhomogeneous metal enrichment. However, this is most 
likely to happen in more massive halos, as smaller halos will 
experience rapid self-enrichment or lose their gas entirely af- 
ter the first Pop III supernovae. Massive halos are also more 
likely to be made up of smaller halos with disparate chemical 
evolution histories (see footnote 9). Such massive halos are 
more common at z ~' 7-9, when the Universe is 0.5 Gyr old. 
This cosmic age is substantially larger than the lOMyr self- 
enrichment timescales of halos discussed above, over which 
the conditions to form Pop III are lost. Additionally, the gain 
in hydrogen-ionizing photon production from Pop III relative 
to Pop II is a factor of order 0.6 up to a few (Tumlinson et al. 
2004). When one considers that this could occur only in larger 
halos for at most lOMyr (at a cosmic age of 0.5-1 Gyr), the 
impact on reionization should be relatively small. We have 
confirmed in separate calculations for such a scenario that this 
is the case. We therefore do not consider simultaneous Pop III 
and Pop n star formation within the same halo in our models. 

Note that in our formalism for using weighted-average ion- 

' The assumption that each new halo forms metal-free stars breaks down at 
z ~ 8-9, when the earliest small halos, presumably metal-enriched from their 
first episodes of Pop III star formation, begin to coalesce into larger objects. 
Most of the scenarios considered here reionize earlier than this; those that do 
not would arguably not follow the standard timeline of chemical enrichment 
anyway, as the bulk of baryons would be tied up in long-lived DSP-phase 
dark stars at redshifts z > 9, leading to some delay in the chemical evolution 
of the Universe. 

'" The Pop III IMF is independent of the presence of dark stars in our 
formalism, and is not designed to reflect the late stages of dark star evolution 
in any way. As noted earlier, dark stars die at the end of the DSNMS phase, 
which is effectively their main sequence. 



izing photon fluxes from Pop III, Pop II and dark stars, at each 
cosmological epoch, we count the time-appropriate Qu at that 
epoch arising from all star-forming halos that have collapsed 
over a range of redshifts starting at z = 20. Reionization is 
defined as the overlap of individual ionized regions of H II, 
when its IGM volume filling factor equals unity (Venkatesan 
et al. 2003). 

5. REIONIZATION HISTORY OF A UNTVERSE CONTAINING DARK 

STARS 

In the NC scenario, where dark stars are fueled by 
gravitationally-contracted dark matter only, we find that the 
reionization history of the Universe is effectively identical for 
all values of the dark star fraction /ds (Table 1). This is de- 
spite the fact that during the first 0.4 Myr of star formation, the 
dark stellar population contributes nothing to reionization, re- 
ducing 2tot to (1 -/ds)GpopIii- In the ensuing DSNMS phase, 
the earlier lack of ionizing photons is compensated for by an 
excess relative to the canonical situation (where /ds = and 
there are no dark stars, only the normal Pop III followed by 
Pop II), because Qdsnms > Qpopiii- Although the net impacts 
of these two effects on the redshift of reionization cancel, one 
might expect a difference in the time-evolution of the ioniza- 
tion fraction relative to the simple Pop 111-1-11 case; due to the 
short durations of the DSP and DSNMS phases in the NC sce- 
nario, this effect is however too small to notice. 

In Fig. 1, we show ionization histories in the extreme cap- 
ture (EC) scenario, for various combinations of /ds and the 
DSP lifetime fosp- Here we see a marked effect on reioniza- 
tion, with populations containing large numbers of dark stars 
and/or relatively long-lived ones resulting in substantially de- 
layed reionization. As expected, the larger the values of /ds 
and fDSP, the larger the effect. For larger /ds and ^dsp, Zrdon 
occurs later, as Qtot remains at lower values for longer periods 
of time during the earliest part of the star-formation history 
in each halo. For longer delays, reionization occurs more 
quickly as IGM ionization becomes able to build up more 
quickly. This is a result of two effects: the main sources of 
ionizing photons turning on at lower redshifts, and the fact 
that the IGM density decreases rapidly with decreasing red- 
shift, leading to increasing IGM recombination timescales at 
later cosmic times. This explains why the rapidity of reion- 
ization becomes more pronounced with increasing fDSP than 
with increasing fos- 

Comparing our results in Table 1 with the constraints from 
WMAP7, we find that a few of our cases can be immedi- 
ately ruled out by simply having too low a value of Tg, the 
optical depth to electron scattering, or Zreion- The global 
fit to WMAP7 and other cosmological data (Komatsu et al. 
2011) implies that the integrated optical depth back to the 
surface of last scattering (i.e. recombination at z ~' 1090), is 
Te ^ 0.088 ±0.014. For a simple step function reionization 
history, this is Zreion = 10.6 ± 1.2. A direct comparison with 
Table 1 would in principle rule out all cases with Tg < 0.074, 
or with Zreion ^ 9.4. Note that the Te constraint is more direct, 
and that limits based on Zieion are only approximate owing to 
the WMAP7 assumption of a simplified reionization history 
(step function ionization at fixed Zieion)- Our models are more 
realistic and track the details of reionization with variations in 
galaxy halo masses and astrophysical parameters, so that our 
derived Zreion and that from WMAP7 may not be directly com- 
parable. Thus, comparing the values of Te from Table 1 to the 
WMAP7 limits, we see that the EC case with fosp = 500 Myr 
and /ds = 1 is ruled out. More EC cases, e.g. fosp = 150Myr 
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Fig. 1. — Ionization histories of a Universe containing dark or semi-dark stellar populations, in the extreme capture (EC) scenario, where dark stars live 
extended lives as cool, diffuse objects. Here we plot the fraction of hydrogen in H II as a function of redshift. Top panels compare histories for different dark star 
fractions /jjs, in situations where dark stars live a long time in the DSP phase (/qsp = 150 Myr and 500 Myr). Bottom panels compare histories for different DSP 
lifetimes, in situations where dark stars make up a substantial fraction of the mass budget of the first population of stars (/ds =0.6 and 1). Longer-lived and more 
numerous EC dark stars delay reionization. Here we have assumed that the fraction of the initial star-forming baryonic mass budget of the Universe that does not 
form dark stars instead forms normal Pop III stars, which die after 10 Myr and are replaced by newly-born Pop II stars. Following the death of the population of 
dark stars, and depending upon their time of death, they are replaced either directly with newborn Pop II stars, or with newborn Pop III stars, which themselves 
later die and get replaced by newborn Pop II stars. 



and /ds = 0.9-1, or fosp = 500 Myr and /ds > 0.75, would 
be ruled out at face value if we were to compare directly with 
the WMAP7 Umit that Zreion = 9.4- 1 1 .8. Similarly, one might 
also conclude that many models with low ^dsp and /ds are 
ruled out for producing z,eion and Tg exceeding the upper limit 
of the WMAP7 error band. Whilst this is indeed true when 
/*/esc = 0.005 as we assume here, such limits are not really ro- 
bust to variations in astrophysical parameters. We discuss in- 
tegrated optical depths and corresponding constraints in more 
detail in the following section. 

For the meager capture (MC) scenario, the effects are less 
dramatic; in Fig. 2 we show a zoomed-in section of the full 
history of reionization, for a few combinations of /ds and the 
DSNMS lifetime fDSNMS- Here, increasing /ds and fDSNMS 
results in progressively earlier reionization, the exact oppo- 



site trend relative to the EC cases. This is because Qdsnms > 
Gpopiii, so extending the DSNMS phase and increasing the 
fraction of baryons contained in it causes reionization to hap- 
pen more quickly than with only a normal Pop 111 IMF. 
Again, as expected the effects increase with increasing /ds 
and fDSNMS ■ The reason the MC scenario has a much smaller 
effect on reionization than the EC scenario is that its dura- 
tion is already much more strongly constrained, in this case 
by the fusion-burning timescale of core hydrogen during the 
DSNMS phase. 

In Fig. 3 we show the impact of varying the product of the 
astrophysical parameters f^, and /esc from our canonical value 
of /t/esc = 0.005 to the extreme values of 0.02 (/t = 0.1, /esc = 
0.2) or 10""* (/* = /esc = 0.01). We give the resulting reioniza- 
tion histories both for a standard Pop III without dark stars. 
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Fig. 2. — Ionization histories of a Universe containing dark stars, in the meager capture (MC) scenario, where dark stars effectively receive a small extension 
to their main-sequence lifetimes. The two panels compare histories for different dark star fractions /ds, where dark stars live in the DSNMS phase for close to 
the maximum time allowed by core hydrogen depletion (/dsnms = 3 and 6 Myr). Longer DSNMS lifetimes and larger dark star fractions lead to small increases 
in the speed of reionization (note the zoomed axes relative to Fig. 1). Non-dark aspects of the calculations are as described for Fig. 1 . 
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Fig. 3 . — Impacts of varying astrophysical parameters upon the reionization history of a Universe containing no dark stars (left), and one containing EC dark 
stars with tjysp = 1 50 Myr, /bs = 1 (right). Here we show the effects of varying the product of the star-forming baryon fraction /* and the ionizing photon escape 
fraction /esc. The variation of astrophysical parameters induces a similar change in the reionization history of the Universe to dark stars (left), but has a slightly 
reduced impact when applied to reionization scenarios that include dark stars (right). Variations in /*/esc can not only delay reionization as EC dark stars do, but 
also speed it up, to a much greater extent than MC dark stars are able to do. 



and for an EC example with /os = l,fDSP = 150 Myr. Here, 
we see that within this range of astrophysical uncertainties, a 
large range of reionization histories is possible. Indeed, in the 
most extreme cases, dark stars have a similar magnitude effect 
as the variation of astrophysical parameters. This degeneracy 
is unfortunate, but not unexpected; substantial uncertainty ex- 
ists in reionization models at present, even before introduc- 
ing the possibility of stellar populations including dark stars. 
Hearteningly however, the impact of the astrophysical uncer- 
tainties is reduced in situations where dark stars play a sig- 
nificant role, as can be seen by comparing the two panels of 
Fig. 3. 



The results we present here agree broadly with those of 
Schleicher et al. (2009), but the correspondence is not imme- 
diately obvious. Schleicher et al. considered reionization from 
'MS-dominated' (main sequence), and 'CD' (capture domi- 
nated) dark stars, roughly corresponding to our own MC and 
EC scenarios, respectively. They investigated the case where 
/ds = 1, showing that the higher ionizing photon fluxes of the 
DSNMS phase hasten reionization, whereas the lower fluxes 
of the DSP phase delay it. This is in good agreement with 
what we show here, and earlier predictions by Yoon et al. 
(2008). Where we differ from Schleicher et al.'s analysis is in 
our atmospheric modeling (we use actual model atmospheres 
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rather than black-body spectra), and in the details of our pop- 
ulation modeling. 

Here we carefully treat the allowed lifetimes of the differ- 
ent dark star phases, taking into account existing limits from 
the timescales of core hydrogen burning and self-annihilation 
of the dark matter halos surrounding dark stars. Schleicher 
et al. (2009) adopted the ionizing photon fluxes of Yoon et al. 
(2008), which accounted for the limit from core hydrogen 
burning, but not from halo self-annihilation or disruption (al- 
though both did at least acknowledge that this might be a 
concern). Many of these lifetimes we now know to be theo- 
retically inaccessible due to self-annihilation constraints (see 
the paragraph discussing viable ranges of fpsp and fosNMS in 
Sec. 3.1). Many of the models Schleicher et al. considered 
were accordingly ruled out by the redshifts of reionization or 
integrated optical depths. Schleicher et al. (2009) also allowed 
only one or the other of the DSP or DSNMS phases in their 
calculations, whereas we include both self-consistently. 

Schleicher et al. (2009) went on to investigate whether the 
effects of dark stars can be compensated for by more com- 
plicated reionization histories, varying assorted astrophysical 
parameters and introducing a second period of recombination 
followed by a late starburst, leading to a two-stage reioniza- 
tion history. Here we have shown that although dark stars 
can have a significant impact upon the reionization history 
of the Universe, they need not necessarily. Even in cases 
with the most extreme effects (e.g. EC scenarios with /ds = 1, 
fosp > 100 Myr), Fig. 3 shows that ad hoc scenarios like those 
considered by Schleicher et al. are not necessary to reconcile 
dark stars with reionization constraints; a simple increase in 
/*/esc does the job quite well enough. 

6. IMPACTS ON THE COSMIC MICROWAVE BACKGROUND 

6.1. Electron-scattering optical depths 

Following ShuU & Venkatesan (2008), for each of our 
reionization histories we calculate the optical depth from the 
present day to a redshift z due to Thomson scattering as 



Te(z) = 



H, 



ne(z)o-T 



Jo 



(i+z) nui+zy+^A' 



-dz, 



(6) 



where 



n.(z)=^(l+z? 
87rGmH 



XfHn(z)+- {/Hen(z) + 2/Hen,(z)} 



(7) 



is the number density of free electrons. Here Hq, il^, iljy 
and r^m are the present-day values of the Hubble constant, 
mass fraction, baryon fraction and dark energy fraction of the 
critical density of the Universe, respectively. The electron- 
photon Thomson-scattering cross-section is given by aj, the 
mass of hydrogen by me, the primordial hydrogen mass frac- 
tion by X, and the primordial helium fraction by F w 1 -X. 
The ionization fractions /hh, /hch and /eem refer to the frac- 
tion by number of hydrogen or helium atoms respectively in 
the ionization states H II, He II and He III. We assume that the 
number of electrons provided by ionization from He I to He II 
directly tracks hydrogen ionization (i.e. /Hen + /Hem = fun), 
leading to 



«e(z) = 



3f7b^() 
SnGmii 



(i+zY 



3Y\ Y 

1-^1 /hii(z)+ jftiemiz) 



(8) 



We assume a simple step-function ionization model for He III, 
with fueiniz > 3) = and /Hem(z < 3) = 1. We also assume 
a residual electron fraction from recombination, present even 
before reionization at the level of Xg = 2.1 x 10"'*. This num- 
ber comes from recombination modeling in CAMB (Sec. 6.2). 
For both the optical depths based on Eq. 6 and CAMB calcula- 
tions, we use the same values for cosmological parameters as 
in our reionization calculations, coming from WMAP 7-year 
results (Larson et al. 201 1). We give optical depths integrated 
up to the surface of last scattering (z ^ 1090) for each of our 
parameter combinations in Table 1 . 

In Fig. 4 we show the evolution of optical depth correspond- 
ing to the ionization histories detailed in Fig. 1. As might be 
expected from the ionization curves, longer-lived and more 
numerous EC dark stars result in smaller electron-scattering 
optical depths, as they reionize the Universe later. The result- 
ing integrated optical depths across the entire EC parameter 
space are summarized in Fig. 5, where we plot Tg as a contin- 
uous function of fosp and /ds ■ 

Similarly, we show a zoomed-in section of the optical depth 
curves for the longest-lived MC case (tqsnms = 6 Myr) in 
Fig. 6. In this case, the smaller variations in reionization 
history have a correspondingly smaller (and opposite) effect 
upon Te, leading to slightly larger optical depths than in the 
/ds = case. 

We also show in Figs. 4 and 6 the la measurement of 
the integrated optical depth to last scattering from WMAP7 
(Te = 0.088 ±0.014; Komatsu et al. 2011), along with a pro- 
jected Planck sensitivity to the same quantity (Colombo et al. 
2009), assuming the two experiments measure the same cen- 
tral value. ' ' Assuming that our chosen astrophysical reioniza- 
tion parameters in our canonical models are correct (/*/esc = 
0.005), a significant part of the more extreme end of the pa- 
rameter space is akeady ruled out at better than a standard 
deviation by the WMAP7 measurement of Tg; much more 
will be excluded by Planck. This is summarized in Fig. 7, 
where we plot the la exclusion curves implied by WMAP7 
and Planck measurements of Tg in the fDSP-/DS plane, for the 
EC scenario. 

Noticeably however, even the standard /ds = 0, Pop IIIh-II 
model exhibits tension with the WMAP7 data at more than 
the la level. As we discuss in Sec. 6.3, this should be taken 
with something of a grain of salt: similar variations in optical 
depth can be produced by reasonable changes in astrophysi- 
cal parameters, so this should not be considered evidence for 
the existence of dark stars at this stage, even at the la level. 
Robust constraints could be obtained by full joint parameter 
scans of cosmological and (dark) reionization models. Whilst 
such an exercise is well beyond the scope of this paper, the re- 
sults we present in this section clearly indicate that dark stars 
can have a significant impact upon the reionization history of 
the Universe, and therefore, the integrated optical depth of the 
CMB. 

6.2. Polarization 

We calculate the effects of the different reionization histo- 
ries on the polarization (EE) CMB power spectra by modi- 
fying the Boltzmann code CAMB'^ (Lewis et al. 2000). In- 

' ' The result we use from Colombo et al. (2009) is in fact the ratio of 
1(7 uncertainties on the measured value of Te obtainable with WMAP5 and 
Planck in the absence of foregrounds. We thus also assume the same percent- 
age degradation in accuracy for both WMAP5 and Planck when mapping 
from expected results without foregrounds to final limits. 

'^ http://camb.info 
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Fig. 4. — Evolution of the lookback optical depth from the present day to redshift z according to Eq. 6, for EC dark stars of varying lifetimes and abundances. 
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Stead of the simple hydrogen reionization model included in 
CAMB, our modified version uses the H ionization fractions 
presented in Sec. 5 as the basis for its reionization calcula- 
tions. We include contributions to the total electron fraction 
from H II, He II and He III with the assumptions stated ear- 
lier, as well as a residual electron fraction. As in the stan- 
dard CAMB reionization calculation, we assume that elec- 
trons from He II track those from H II, and model contribu- 
tions from He III with a smoothed step function centered at 
z^ 3.5. We take the residual electron fraction after recom- 
bination to be Xe = 2.1 X 10""*, based on the output of REC- 
FAST within CAMB. Using CAMB's highest accuracy set- 
ting, we calculated the temperature (TT), polarization (EE) 
and cross (TE) power spectra, as well as integrated optical 
depths. We checked that the optical depths computed with 
CAMB agree to within their stated numerical accuracy with 
those from Eq. 6 (as presented in Table 1 ), and verified that the 



slight difference in the treatment of He III here and in Sec. 6. 1 
has a negligible impact upon integrated optical depths. 

In Fig. 8 we plot EE polarization spectra for the same EC 
cases as illustrated in Figs. 1 and 4. The spectra are normal- 
ized to the corresponding EE spectrum obtained in the stan- 
dard /ds = 0, Pop 111-1-11 scenario, in order to investigate the 
ability of CMB experiments to distinguish dark stars from 
standard reionization. To this end, we also plot the uncer- 
tainty on the normalization due to cosmic variance, and the 
combination of cosmic variance and total readout noise ex- 
pected across the 70, 100 and 143 GHz channels in the first 
14 months of Planck operation (Colombo et al. 2009). Large 
parts of the parameter space are distinguishable from /ds = 
in a cosmic -variance-limited experiment, and a number of the 
more extreme scenarios are even detectable by Planck at bet- 
ter than \g. Although essentially all such models may be dis- 
favored anyway by Planck's measurement of the integrated 
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optical depth, the EE spectrum would nonetheless provide a 
complementary (albeit weak) statistical handle via which to 
increase the power of full parameter scans to exclude such 
dark star models. 

We do not show TT or TE spectra for the EC scenario, as 
they exhibit less striking deviations from the corresponding 
spectra of the standard Pop lU-i-II scenario at low multipoles / 
(large angular scales) than the EE curves do. We do point out 
however that the TT and TE spectra exhibit damping at large 
I due to the changing optical depth, which is more clearly 




DSP lifetime tosp (Myr) 

Fig. 7. — Implied \a detection/exclusion regions in the /ds^'dsp plane 
for EC dark stars, based on the integrated optical depth to last scattering ob- 
served by WMAP7 (Komatsu et al. 2011). We also show projected Planck 
constraints (Colombo et al. 2009), assuming that the same central value 
(Te = 0.088) is measured by Planck as by WMAP7. To guide the eye, the 
depth of shading is proportional to the likelihood of each parameter com- 
bination. For a product /*/csc = 0.005 of the star-formation efficiency and 
UV photon escape fraction, parameter combinations outside the red WMAP 
shaded region are excluded at greater than 1 cr by existing data. Combinations 
outside the shaded blue region will be excludable at better than 1 cr by Planck. 
Note however that variations in /*/esc will shift these regions substantially; 
refer to discussions in the final paragraph of Sec. 6.1 and in Sec. 6.3. The 
same interpolation methods were employed in this figure as in Fig. 5; the 
slight upturn of the boundaries of the Planck region at large fusp is again an 
artifact of the interpolation. 

visible than in the EE spectra. We also do not show power 
spectra for the MC or NC cases, as they show little deviation 
in general from the standard Pop lll-i-ll case. 

6.3. Astrophysical uncertainties and implications for 
parameters of reionization models 

In Fig. 9 we show the impacts of varying astrophysical pa- 
rameters upon CMB observables. Here we give the variations 
in optical depth and EE polarization resulting from the same 
variations of /*/esc as in Fig. 3. For EE spectra, we nor- 
malize to the corresponding curve with default astrophysical 
parameters in each case; i.e. to the standard /t/esc = 0.005 
Pop IIIh-II spectrum for the /ds = curves, and to the /os = 1, 
fosp = 150 Myr, /*/esc = 0.005 spectrum for the curves where 
/ds = 1, fosp = 150 Myr Comparing with Figs. 4 and 8, we see 
again that variations in the astrophysical parameters within 
reasonable ranges can have similar effects (both in strength 
and character) to variations in the dark star parameters. Al- 
though this means that the impact of dark stars on the CMB is 
very difficult to unambiguously disentangle from existing the- 
oretical uncertainties in reionization modeling, it also serves 
as a clear indication that the potential effect of dark stars upon 
CMB observables affected by reionization could be quite sig- 
nificant. 

In some cases, it is possible that specific regions of the 
reionization parameter space that are ruled out in standard 
Pop 111-i-ll-only scenarios by the current WMAP7 (or pro- 
jected Planck) data, are reopened by the possibility of having 
dark stars. For example. Fig. 9 reveals that the extreme cases 
of varying astrophysical parameters (/^/esc = 0.02 or 10"'*) are 
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ruled out for the Pop Ill+II-only scenarios by current CMB 
data. However, the case of /*/esc = 0.02 is not ruled out for 
the EC dark star case with /ds = 1, fosp = 150 Myr. Having 
more dark stars in the EC scenarios mimics a reduction in 
the astrophysical efficiency - i.e. increasing /ds allows for 
greater astrophysical efficiency (larger /t, larger /esc or both) 
in reionization models. This is potentially an important fac- 
tor to consider when placing constraints on the astrophysical 
aspects of reionization from CMB data (Haiman & Holder 
2003; Shull & Venkatesan 2008). 

The impacts of uncertainties in cosmological parameters on 
Te are subdominant to those from astrophysical uncertainties. 
For the same two scenarios shown in Fig. 9 (no dark stars and 
EC dark stars with /ds = 1 and fosp = 150 Myr), we find that 
varying a% over the 1 a allowed range of the WMAP7 fit we 
used as input for all calculations (cts = 0.801 ±0.030), leads 



to a change in the redshift of reionization of less than 1, and 
a change of no more than 0.01 in t^. This is relatively small 
compared with the variations in Fig. 9. 

7. CONCLUSIONS 

We have calculated reionization histories, CMB optical 
depths and anisotropy power spectra for a broad range of stel- 
lar population models containing dark stars. We identified 
three distinct regimes: no capture (NC), where dark stars live 
as cool, semi-diffuse objects for no more than ^ 0.4 Myr and 
then quickly move on to the main sequence; meager capture 
(MC), where dark stars undergo the same initial evolution as 
in the NC case, but exhibit slightly decreased surface tem- 
peratures and moderately increased lifetimes on the main se- 
quence, and extreme capture (EC), where dark stars live an 
extended life as cool, semi-diffuse objects before reaching the 
main sequence. 
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NC dark stars have effectively no impact on reionization or 
its signatures in the CMB. MC dark stars cause the Universe 
to reionize more quickly the longer-lived and more numerous 
they are, producing a slight increase in the redshift of reion- 
ization. This leads to small increases in the total integrated 
optical depth of the CMB. 

EC dark stars can have dramatic and completely opposite 
effects to MC dark stars, delaying reionization to as late as 
z ^ 6. Greater numbers of EC dark stars and longer dark star 
lifetimes lead to increased delays, and correspondingly de- 
creased integrated optical depths to the last scattering surface 
of the CMB. Using WMAP7 observations of the integrated 
optical depth, we have been able to rule out more extreme ar- 
eas of the parameter space covered by the EC scenario. Planck 
will improve these bounds significantly. EC dark stars also 
produce a characteristic suppression of EE polarization power 
on the largest scales of the CMB, as well as a slight enhance- 
ment at small scales due to the decreased optical depth. 

Many of these effects can be mimicked or compensated for 
by changes in the astrophysical parameters of standard reion- 
ization models. Disentangling the impact of dark stars from 
other theoretical uncertainties in reionization modeling will 
be challenging, even with Planck. Not only can dark stars 
have a substantial impact on reionization and its signatures 
in the CMB, but the addition of dark stars to standard reion- 
ization models can in fact substantially increase the range of 



astrophysical parameters that can be made consistent with ex- 
isting (and future) observations. 
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